extern doublereal snrm2_(integer *n, real *x, integer *incx);
